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Overview 

Problem:  The  mixing  of  differently-sized  particles  is  already  difficult  due  to  such  phenomena  as  the  Brazil  nut 
effect.  Cohesion  between  particles  (e.g.  arising  from  van  der  Waals,  capillary,  or  electrostatic  forces)  adds  a  level 
of  complexity  to  the  mixing. 

Agglomerates  need  to  be  broken  up  to  achieve  good  mixing 

Mixing  harder  to  break  up  agglomerates  can  lead  to  greater  and  faster  “unmixing”  (e.g.  the  Brazil  nut 
effect) 

While  the  strength  of  cohesion  can  often  be  altered  (e.g.  via  surfactants),  the  effects  of  doing  so  on  the 
mixing  process  have  not  been  fully  explored  and  quantified,  especially  when  the  particles  feature  a  large  size 
range 

To  control  and  maximize  the  homogeneity  of  mixtures  of  cohesive  particles,  we  seek  to  develop  a  model  that 
relates  a  mixture’s  homogeneity  to  its  components’  properties  and  mixing  conditions 


Approach 

•  Simulate  the  mixing  process  of  cohesive  particles  with  a  variety  of 
system  parameters  to  determine  their  effect  on  mixtures’ 
homogeneity 

Simulation  technique:  DEM  (Discrete  Element  Model/Method) 

Replace  “complex  mixing  procedures”  with  simple  plane  shear  flow 

Examine  only  one  species,  i.e.  all  particles  have  same  density,  Young’s 
modulus,  roughness,  etc.  -  only  size-segregation  and  agglomeration  due 
to  cohesion  will  be  examined 

Simplify  size  and  shape  distribution  by  using  spherical  particles  of  2  sizes 
with  7: 1  diameter  ratio 

*  Individual  small  particles  are  able  to  fit  in  between  packed  big  particles 

*  Homogeneity  can  have  a  large  effect  on  total  packing  fraction 

Properties  to  examine 

*  Particle  cohesiveness  ,  or  “stickiness”,  which  can  be  controlled  in  real  life  by 
surfactants  -  some  or  all  particles  can  be  altered 

*  Shear  rate 

*  Interstitial  fluid  (e.g.  binder)  -  buoyancy  (fluid  density)  and  viscosity  (not 
discussed  here) 
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DEM  simulations* 

The  particular  DEM  code  is  a  modification  of  LAMMPS,  an  open 
source  code  from  Sandia  National  Labs 


•  Every  particle  simulated  as  a  discrete  spherical  volume  w /  mass 

•  Particles’  positions  and  velocities  calculated  by  integration  of 
Newton’s  2nd  law  (F  =  ma) 


•  Particles  are  allowed  to  overlap  (“soft”)  and  impose  forces 
analogous  to  inelastic,  frictional  springs  when  they  collide 


Fta„  =  min(kjvtdt,  nFnormal) 


tan 

tangential 

spring 

constant 


I 


friction  coefficient  — 
limits  tangential  force 
according  to  Coulomb’s 
law  of  friction 


•  Cohesion  between  particles  can  be  simulated  by  allowing  them  to 
attract  each  other  via  the  van  der  Waals  force 

A  is  varied  for  each 

Fvan  der  Waals  =  -Aefl<W24s2  gr°UP  of  particles  to 

,  ..  investigate  influence  of 

Hamaker  constant  surface  separation  ° 

cohesion/surfactants 


*  P.  A.  Cundall  &  D.  L.  Strack, 
Geotechnique  29  (1979) 
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Simulation  details 


45679  small  (red)  particles  of  diameter  d  (~l/4  of  solid  volume) 
367  large  (blue)  particles  of  diameter  7d  (~3/4  of  solid  volume) 

Big  and  small  particles  are  initially  separated 


Ceiling  can  — 
move  vertically 
to  allow  for 
bed 

contraction 
and  expansion 
s.t.  it  applies 
constant 
pressure  Pd/k 
=  2.87  x  10  4 


Systems  bounded  on  top  and  bottom  with  walls  made  of 
small  particles  which  are  moved  to  generate  shear 
Streamwise  &  neutral  directions  are  periodic  (49d  in  length) 

Interstitial  fluid  (binder)  is  not  explicitly  simulated 
It  is  assumed  to  fill  in  all  empty  spaces  and  only  exerts  a  buoyant 
(upward)  force  on  all  particles  (except  the  bounds):  Fb  =  g/^V 

H  5 


Investigating  the  influence  of  the  binder  density 


Case  I :  Particles  neutrally  buoyant  in  binder 


Pb  Pp 


°  Buoyant  forces  negate  gravitational  forces  and  the  resulting  gravity- 
driven  segregation  (e.g.  Brazil  nut  effect) 

°  Properties  “after  long  times”  (steady-state?)  are  of  interest 


°  If  there  is  no  cohesion,  then  there  will  be  no  driving  force  for 
segregation  ->  “best  case” 

*  As  cohesion  of  all  particles  is  increased,  segregation  should  increase 

*  Increasing  shear  rate  should  lead  to  agglomerates  breaking  up  and  small 
particles  getting  into  gaps  between  big  particles,  improving  homogeneity 


*  The  effects  of  altering  the  cohesiveness  of  just  one  group  of  particles  is  not 
immediately  apparent 


Case  2:  Particles  denser  than  binder 


l-8pb  =  Pp 


°  Big  particles  will  rise  to  top  over  time 

*  How  homogeneous  do  mixtures  get  before  this  happpens? 

°  Unknown  if  and  how  varying  shear  rate  and  particles’ 
cohesiveness  affects  homogeneity  throughout  mixing 

°  Initial  conditions  may  have  strong  influence  on  mixing 
quality  -  particles’  paths  will  be  affected  by  gravity 
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Simulation  dimensionless  parameters 


When  pb  =  pp,  particles  in  effect  do  not  feel  any  gravitational  force  (g  is 
unimportant),  so  g  is  not  used  to  scale  cohesiveness  or  shear  rate 


P  vdW 

k  =  normal  stiffness 
of  small  particle 


Atdi 


24  kds;„. , 


=  2X10  ~5d 


“modified  Bond  number"  - 
dimensionless  particle  cohesiveness 


minimum  surface  separation  used  to  prevent  van  der 
Waals  force  model  from  diverging  during  collisions 
Assuming  d  »  20pm,  smin  «0.4nm,  corresponding  to  typical 
intermolecular  distances 


To  eliminate  the  particle  size  dependence  of  Bo*  and  to  rescale  it  to  be 
on  average  of  order  unity,  we  define  a  scaled  modified  Bond  number: 


Bo**  =  4000 Bo* 


S  =  Bo 


** 

small 


B  -  Bo 


** 

big 


Shear  rate  expressed  using  “scaled  stiffness”  - 
dimensionless  inverse  square  shear  rate 


k*  = 


k 


3-2 


Pbdf 


t 

Binder  density  is  kept  the 
same  among  all  simulations, 
unlike  the  particle  density 
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Order  metrics 


•  Spatial  variance  in  volume  fraction  for  each  particle  size 

°  “Most  straightforward”  measurement  of  homogeneity 

•  Total  volume  fraction 

°  If  well-mixed,  small  particles  fill  in  gaps  between  large  particles,  leading 
to  a  high  volume  fraction 

°  However,  a  low  volume  fraction  does  not  necessarily  mean  poor  mixing 

Estimating  the  average  size  of  clusters  of  small  particles-  smaller 
cluster  size  means  better  mixing 

°  Model  has  recently  been  proposed  to  estimate  average  cluster  size, 
essentially  based  on  radial  distribution  function 


d 


cluster 


=  2 r(L{rag)  +  rag) 


S.  Gallier,  Propellants  Explos. 
Pyrotech.  34  (2009) 


Besag  L-function 

J.  Besag,  J.  R.  Stat.  Soc.  B39  (1977) 


Maximum  value  of  L  [r ) 


occurs  at  r  =  rag  (aggregate  radius) 


Ripley  K-function  (g(r)  integrated  over  space) 
B.D  Ripley, J.Appl.  Probab.  13  (1976) 
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Case  I  Results: 
Particles  neutrally 
buoyant  in  binder 

Pb  “  Pp 


***  The  results  presented  in  this  section  are  from  simulations  that  have 
been  run  for  “long  times,”  not  necessarily  to  steady-state  *** 


As  expected,  when  particles  are  non-cohesive, 
the  mixtures  appear  homogeneous 


k*  =  109  k*=  108  k*=  107 
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When  small  particles  are  not  very  cohesive  (S  =  0.1)  , 
the  composites  also  appear  relatively  homogeneous 


e  =  o.oi 


k*=  107 


k*=  io8 


6  =  0.1 


Small  particles  do  not 
agglomerate  enough  to 
prevent  them  from  mixing 


8=1 


When  small  particles  are  more  cohesive  (S  =  l),the  composites 
appear  most  inhomogeneous  at  low  shear  rates  and  at  very  high 
and  very  low  big-particle  cohesion 

8  =  0.1  8  =  1  8  =  10 


k*  =  107 


k*  =  108 


Big  particles 
pack  more  tightly 
-  small-particle 
agglomerates 
cannot  fit  in  gaps 


k*  = 


Increased 
cohesion  inhibits 
small  particles’ 
ability  to  enter 
gaps  between  big 
particles 
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normalized  st.  dev.  in  small  ptl.  vol.  fract. 


Spatial  variance  in  small-particle  volume  fraction 
more  or  less  reflects  what  was  seen  in  the  images 


normalized  st.  dev.  in  big  ptl.  vol.  tract. 


Spatial  variance  in  big-particle  volume  fraction  still 
reflects  images  but  behaves  somewhat  differently 

0.4 
0.35 
0.3 
0.25 
0.2 
0.15 
0.1 
0.05 

10'  10“  10 


1  1  — 1 — 

★  no  cohesion 

S  =  0.1,  B  =  0.01 
-  A  S  =  0.1,  B  =  0.1 
■  S  =  0.1,  B  =  1 

#  S=  1,  B  =  0.1 

S=  1,  B=  1 

V  S=  1,  B=  10 

- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

A 

Order  ^ 

*  reversed 

Significantly 
y  larger 
/  difference 

- 

f 

Inhomogeneous? 

Homogeneous? 

▼ 

■  ★ 

l 

_J _ 1 _ 1 _ 

*  ■ 

_ i _ i _ i _ i _ i _ i _ 1 _ i _ i _ i _ i _ i _ i _ i _ i _ 1 _ 
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cluster  size  (d) 


Cluster  size  measurements  also  reflect  mixing  quality 
observed  in  images  but  again  manifest  different  behavior 

4 

3.5 

3 

2.5 

2 

1  5 

„  rJ  „  -8  .  -9 

10  10  10 


^  i  i 


i  r - i  !  T - 1  i 


★  no  cohesion 

S  =  0.1,  B  =  0.01 
AS  =  0.1 ,  B  =  0.1 
■  S  =  0.1,  B  =  1 

♦  S  =  1,  B  =  0.1 
S  =  1,  B  =  1 

▼  S  =  1,  B  =  10 


v 


Why  are  data 
separated  by 
so  little? 


♦ 


i 


Jnhomogeneous? 

Homogeneous? 


_l _ I _ I _ I _ l_ 


_l _ I _ I _ I _ I _ I _ I _ 
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solid  volume  fraction 


Volume  fraction  does  not  capture  the  homogeneity  well 
but  is  nonetheless  greatly  affected  by  particle  cohesion 

0.8 

0.75 

0.7 

0.65 

10  10  10 


— I - 1— 


— I - 1 - 1 - 1 - 1 - 1— 


* 


★  no  cohesion 

iS  =  0.1,  B  =  0.01 
S  =  0.1,  B  =  0.1 

JS  =  0.1,  B  =  1 
S  =  1,  B  =  0.1 
S  =  1,  B  =  1 
▼  S  =  1,  B  =  10 


★ 


★ 

* 


f 


Volume  fraction  decreases 
with  shear  rate  and 
cohesion,  more  strongly 
with  S  than  with  B 


_l _ I _ I _ I _ I _ l_ 
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Steady  state  seems  elusive  in  most  cases 
Should  we  expect  and  wait  for  steady  state? 

Steady  state  has  only  been  reached  under  conditions  that  have 
produced  microstructures  that  look  homogeneous,  but  not  all 
mixtures  that  appear  homogeneous  are  at  a  steady  state 

If  we  continue  to  run  these  simulations,  will  the  inhomogeneous 
mixtures  eventually  become  homogeneous  or  will  they  stay 
inhomogeneous? 

°  If  they  do  eventually  become  homogeneous,  simulations  would  take 
unfeasibly  long  to  reach  homogeneous  states  at  current  pace 

°  Should  the  simulations  simply  be  stopped,  and  if  so,  when? 

The  spatial  variances  in  the  concentrations  of  the  different 
particles  usually  continue  to  evolve  with  time  well  after  the  total 
volume  fraction  and  average  small-particle  cluster  size  have 
reached  steady  state 


Case  2  Results: 
Particles  denser 
than  binder 

l-8pb  =  Pp 
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When  initially  horizontally  segregated,  large  particles 
generally  reach  the  top  before  good  mixing  can  be 
achieved,  regardless  of  cohesion  and  shear  rate 


k* =  1 07  k*=W  B~_^  8  =  0.1  8  =  0.1 

**-10  k*  =  1 07  k*='\09 


S  =  0.1  S=  1  S  =  1  S=  1  S=  1 

8=1  8  =  1  8=1  8=10  8=10 

k*=  109  k*=W  k*=  109  /c*=  107  k*=  109 


So  we  turn  our  attention  to  simulations 
starting  with  vertically  segregated  particles 
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volume  fraction 


Volume  fraction  reaches  a  peak  when  images  show 
well-mixed  and  drops  off  as  un-mixing  occurs 


s  =  o.i 
8=1 
k*=  107 


0.72  i - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


0  54  I _ L _ _ L _ ' _ L_ _ 1 _ 1 _ 1 _ ,! _ 

0  20  40  60  30  100  120  140  160  130  200 

strain 
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standard  deviation  in  volume  fraction 


Minima  in  the  spatial  variance  in  species  concentration 
may  also  be  used  to  identify  and  quantify  well-mixed  state 


1.5 


S  =  0.1 
B=  1 
k*=  107 
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0.5 
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100 
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standard  deviation  in<)> 


However,  minima  in  the  volume  fraction  standard 
deviation  may  actually  be  achieved  during  un-mixing 


big  particles 
small  particles 
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50 


100 


strain 
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well-defined 


How  useful  is  this  statistic  to  identify 
and  quantify  homogeneous  mixing? 


150 
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The  use  of  the  average  small-particle  cluster  size  to 
identify  and  quantify  homogeneous  mixing  is  questionable 
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Minima  are  generally 
not  well-defined  and 
do  not  necessarily 
correspond  to  states 
of  good  mixing 
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Initial  conditions  may  affect  when  the  most 
homogeneous  mixtures  are  achieved,  but  not  necessarily 
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Not  enough  simulations  have  been  run  at  this  time  to 
describe  a  relation  between  homogeneity  and  cohesion 
and  shear  rate,  but  the  results  obtained  so  far  suggest 
that  it  is  complex  and  may  be  impossible  to  predict 


Summary 

•  DEM  simulations  of  bidisperse  cohesive  particles  in  the  presence  of 
an  inviscid  binder  under  shear  were  performed  to  determine  how 
the  homogeneity  of  the  materials  vary  with  the  particle 
cohesiveness,  mixing  speed,  and  binder  density 

When  particles  are  neutrally  buoyant  in  binder 

Homogeneous  mixtures  are  achieved  when  the  shear  rate  is  large  or  the 
small  particles  are  not  very  cohesive 

When  the  mixtures  are  inhomogeneous,  they  are  less  so  when  the  big 
particles  are  as  cohesive  as  the  small  ones 

The  volume  fraction  fails  to  capture  the  homogeneity  as  well  as  the 
other  metrics  which  all  display  different  behaviors 

When  particles  are  denser  than  the  binder  (and  particles  are 
initially  segregated  vertically) 

A  spike  in  the  volume  fraction  occurs  when  “good  mixing”  occurs, 
followed  by  un-mixing 

The  average  cluster  size  and  spatial  variance  in  particle  species 
concentration  are  not  necessarily  minimized  when  mixtures  are  visually 
homogeneous  and  so  may  not  be  the  best  way  to  quantify  homogeneity 

Initial  conditions  have  a  greater  influence  on  the  temporal  evolution  of 
the  mixtures  than  on  how  homogeneous  the  mixtures  become 
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